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PREFACE 


Photothermal  ionization  spectroscopy  is  an  extremely  sensitive 
technique  for  the  detection  and  identification  of  small  concentrations 
of  impurities  in  semiconductors.  It  has  been  used  extensively  in  Si,  Ge 
and  GaAs.  Identification  of  many  of  the  shallow  donors  in  GaAs  with 
their  corresponding  spectral  peaks  has  been  accomplished.  This  provides  a 
useful  tool  for  the  characterization  of  high  purity  semiconductors  giving 
the  crystal  grower  information  about  type  and  relative  concentration  of 
residual  impurities  in  high  purity  material.  We  describe  here  an 
apparatus  which  has  been  constructed  to  perform  such  measurements  and 
the  theory  of  its  operation.  Spectra  of  samples  from  a  number  of  university, 
industrial  and  government  laboratories  have  been  measured.  These  samples 
have  been  grown  by  a  wide  variety  of  growth  techniques  including  MO-CVD, 

MBE,  LPE  and  both  hydride  and  chloride  VPE.  The  spectra  of  these  samples 
are  shown  in  the  figures  and  analyzed  in  detail  in  the  text.  Finally,  a 
computer  investigation  of  the  effects  of  the  Fourier  transform  of  the 
data  on  spectral  line  shape  is  described.  This  is  of  interest  because 
the  spectral  line  shapes  may  be  useful  as  a  tool  to  study  impurity- 
impurity  interactions.  It  is  thus  imperative  that  the  Fourier  transform 
of  the  data  reflect  the  true  photothermal  ionization  line  shape. 
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BASIC  THEORY 


For  an  isolated  impurity,  a  simple  hydrogenic  model  can  be  used  to 
estimate  the  impurity  energy  levels.  A  donor  atom,  relative  to  the  host 
atom  it  replaces,  has  one  extra  valence  electron  and  one  extra  positive 
charge  on  its  ion  core.  To  the  extent  that  electrons  in  the  conduction 
band  behave  as  free  electrons,  this  is  analogous  to  a  hydrogen  atom  in  a 
vacuum.  Using  the  dielectric  constant  and  effective  mass  appropriate  for 
the  semiconductor  (instead  of  for  the  vacuum)  a  model  Hamiltonian  for  the 
donor  electron  can  be  written: 


(1) 


e  r 
o 


*  2  e 

Except  for  the  replacement  of  m  by  m  and  e  by  — ,  this  is  the  familiar 

o 

Hamiltonian  of  the  hydrogen  atom.  The  replacements  simply  rescale  the 
energy  and  distance  scales  of  the  hydrogenic  solutions  giving  the  energy 
spectrum  and  Bohr  radius  as: 


*4  222  *  22 

(2)  E  =  m  e  /2e  R  n  =  [(m  /m)/e  n  ]  13.6  eV 

no  o 

(3)  ag  =  (R2/m  e2)eQ  =  [eo/(m  /ra) ]  .529  A 


where  (m  /ra)  is  the  effective  mass  ratio  and  e  is  the  static  dielectric 

o 

constant  of  the  semiconductor. 

The  ionization  energies  and  Bohr  radii  calculated  from  (2)  and  (3) 
for  donors  in  several  semiconductors  are  given  in  Table  I.  This  picture 
for  electrons  can  be  more  rigorously  derived  by  approximations  to  one  electron 
band  theory''-  and  can  be  extended  to  the  cases  of  anisotropic  effective  masses 
and  multiple  conduction  band  minima.  When  this  is  done, the  analogue  of  (1) 
becomes  intractable  analytically,  but  the  variationally  derived  energies 
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TABLE  I 


Ionization  Energies  and  Bohr  Radius 
Calculated  from  the  Simple  Hydrogenic  Model 


Semiconductor 

m_* 

m 

eo 

ed 

(meV) 

aB 

(£) 

aB 

a 

0 

a 

(£) 

Si 

0. 4a 

11.40 

41.7 

15 

2.8 

5. 43072 

Ge 

0. 2a 

15.36 

11.5 

41 

7.2 

5. 65754 

GaAs 

0. 0665 

12.5 

5.77 

99 

17.6 

5. 65315 

InSb 

0.013 

16 

.689 

650 

100 

6.47877 

Is— 3p(colc.)  >  - 

E 

g 

- Is— 4p(colc.)  « 

- Is— 5p(calc.) 

I - Is— Conduction  Band  (calc.) 


T  *  4.2*K 
Nq>  48XIO,3cm 

N  »2.IXlO,3cm 


_ JJ _ JSHP? 


agree  well  with  those  deduced  from  absorption  or  photothermal  ionization 
measurements,  at  least  for  the  excited  states. 

Another  correction  to  the  hydrogenic  model  is  due  to  the  fact  that 
the  conduction  band  is  not  precisely  parabolic  about  its  minima.  The 
non-parabolic  corrections  to  the  donor  energy  levels  can  be  estimated 
from  Kanes  k*p  calculation  of  E(k)  for  small  k  and  are  less  than  .005  E^ 
in  GaAs.  The  agreement  between  transition  energies  predicted  from  the 
simple  hydrogenic  theory  and  those  observed  in  a  high  purity  GaAs  sample's 
photoconductivity  spectrum  can  be  seen  in  Fig.l. 

It  is  important  to  note  that  the  hydrogenic  model  above  predicts 

identical  ionization  energies  for  different  donor  species.  This  is  not 

observed  experimentally.  The  reason  for  this  is  that  for  distances  less 

than  a  few  lattice  constants  the  potential  felt  by  the  donor  electron  is 

no  longer  the  simple  screened  caulombic  potential  in  (1)  but  rather  is 

determined  by  the  species  dependent  details  of  the  electronic  structure  of 

the  donor  ion  core  and  the  local  distortion  of  the  lattice  around  it.  To 

our  knowledge, no  satisfactory  calculations  of  this  effect  have  been  made. 

In  spite  of  this  interesting  complication,  for  small  binding  energies  the 

donor  wavefunctions  are  quite  extended  and  their  overlap  with  the  region 

where  the  potential  deviates  from  coulombic  is  small.  The  result  is  that 

the  difference  between  the  predicted  hydrogenic  ground  state  energy  and 

the  measured  value  (termed  the  central  cell  correction  or  chemical  shift) 

is  small.  In  GaAs  the  donor  states  are  sufficiently  extended  that  the 

central  cell  corrections  are  only  a  few  percent  of  the  donor  ground  state 

energy.  The  impurity  species  dependence  of  the  central  cell  corrections 

is  of  great  interest  since  it  has  allowed  the  tentative  identification  of 

3 

many  of  the  donors  in  GaAs.  Measured  ground  state  central  cell  corrections 
for  various  donor  species  in  GaAs  are  given  in  Table  II. 
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*  - 1 

TABLE  II  Central  Cell  Corrections  to  the  Parabolic  Band  Ry(  =  46.  10  cm  ), 
for  Four  Residual  Donors  in  High  Purity  GaAs 


Donor 

Central  Cell  Correction 

Ionization  Energy3’*3 

1 

.064  meV 

5.  800  meV 

(.52  cm'1) 

(46.79  cm'1) 

2 

.  081  meV 

5.817  meV 

(.  65  cm*1) 

(46.  92  cm'1) 

3 

.117  meV 

5.  854  meV 

(.94  cm'1) 

(47.22  cm'1) 

4 

.  200  meV 

5.937  meV 

(1.61  cm'1) 

(47.89  cm'1) 

a.  The  ionization  energies  include  the  deepening  of  the  Is  state  due  to  both 
central  cell  shifts  and  conduction  band  nonparabolicity. 

b.  The  magnitudes  of  the  ionization  energies  are  accurate  within  about  ±  0.  15  cm' 1 
while  the  differences  between  the  energies  are  known  to  about  ±  0. 01  cm'1. 
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EXPERIMENTAL  METHOD 


A  means  by  which  the  energy  level  spectra  and  the  associated  central 
cell  corrections  may  be  measured  in  high  purity  semiconductors  is  photo- 
thermal  ionization  spectroscopy.  This  an  extremely  sensitive  method  for 
analysis  of  impurities  in  high  purity  semiconductors.  Other  methods,  such 
as  optical  absorption  and  mass  spectrometry,  are  inadequate  to  identify 
impurities  at  the  extremely  small  concentrations  present  in  high  purity 
material. 

The  photothermal  ionization  technique  is  based  on  the  fact  that  the 
extrinsic  photoconductivity  spectrum  under  certain  conditions  is  a  line 
spectrum  which  corresponds  to  the  bound  electron  energy  levels  of  the 
impurities.  The  technique's  sensitivity  is  due  to  the  fact  that  the 
magnitude  of  the  photoresponse  does  not  decrease  with  decreasing  impurity 
concentration  down  to  very  low  concentrations.  This  is  a  consequence  of 
the  voltage  sensitivity  of  the  extrinsic  photoconductivity  being  proportional 
to  An/n,  where  n  is  the  total  carrier  concentration  without  illumination 
and  An  is  the  change  in  carrier  concentration  under  illumination.  Both  n 
and  An  are  proportional  to  the  concentration  of  major  and  compensating 

4 

impurities  and  so  this  dependence  cancels  out  of  the  voltage  sensitivity. 

This  feature  of  the  method  enables  the  detection  and  identification  of 

impurity  species  even  in  the  highest  purity  semiconductor  material  which 

,  5 

can  be  grown. 

Photothermal  ionization  spectroscopy  requires  a  tunable  far  infrared 
light  source  and  some  means  for  recording  the  sample's  electrical  conduc¬ 
tivity.  The  sample  temperature  is  chosen  such  that  kT<<E^  and  nearly  all 
of  the  uncompensated  donors  are  occupied  by  a  ground  state  electron.  When 


Che  far  infrared  light  source  is  turned  on,  if  the  photon  energy  hv>E^,  some 
of  these  electrons  can  be  excited  directly  to  the  conduction  band  increasing 
the  conductivity  of  the  sample.  (See  Fig.  2(a).)  In  principle,  it  would  be 
possible  to  measure  ionization  energies  by  measuring  the  threshold  of  such  a 
process,  but  in  practice  the  threshold  is  not  sharp  enough  to  permit  this 
measurement . 

Under  certain  conditions,  however,  the  extrinsic  photoconductivity  is 

dominated  by  a  multistep  process  called  photothermal  ionization  and  it  is 

this  process  through  which  the  impurity  energy  spectra  are  obtained.  The 

mechanism  for  photothermal  ionization  can  be  understood  by  considering  a 

donor  electron  in  its  ground  state.  This  electron  can  make  a  transition  to 

one  of  the  donor's  excited  states  by  absorption  of  a  far  infrared  photon 

of  appropriate  energy.  The  still  bound  electron  can  then  be  excited  to  the 

conduction  band  (by  absorption  of  phonon(s))  where  it  can  contribute  to  the 

samples  conductivity.  (See  Fig.  2(b).)  Suppose  one  is  probing  the  ls-2p 

transition.  When  the  light  source  is  tuned  through  hv  =  (E2p-Els^  there 

is  a  peak  in  the  electrical  conductivity  which  corresponds  to  a  ls-2p 

photon  absorption  followed  by  absorption  of  phonon(s)  which  excites  the 

impurity  electron  into  the  conduction  band.  The  sample  temperature  is  chosen 

such  that  there  are  ample  phonons  energetic  enough  to  carry  an  electron  from 

the  excited  state  of  interest  to  the  conduction  band,  and  yet  few  phonons 

sufficiently  energetic  to  excite  electrons  from  the  impurity  ground  state. 

£ 

In  GaAs  this  limits  the  temperature  to  a  range  of  about  1.5  to  6  Kelvins. 

When  a  GaAs  sample  containing  shallow  donor  impurities  is  placed  in  a 
magnetic  field,  the  three-fold  degeneracy  of  the  p  states  is  lifted.  These 
states  are  labeled  in  the  same  way  as  the  corresponding  states  in  the  hydrogen 
atom  (np(m*0,±l) ) .  When  the  magnetic  field  is  increased  from  0  to  7.5  kilogauss. 
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Figure  2. 


Figure  2  (b)  \ 

I 

Microscopic  picture  of  photoconductive  processes. 
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the  dominant  ls-2p  peak  seen  in  the  zero  field  spectrum  splits  into  three 
peaks  corresponding  to  ls-2p  (m=-l,0,+l)  transitions.  (See  Figs.  3(a)  and 
(b).)  A  triplet  also  appears  at  higher  energy  which  is  attributed  to  the 
ls-3p(m=  -1,  0,  +  1)  transitions. 

As  the  field  is  increased  further, the  lines  in  the  photoconductivity 
spectrum  continue  to  separate  and  become  narrower.  (See  Figs.  3(c)  and  (d).) 
The  width  and  shape  of  these  lines  are  determined  by  Stark  shifts  of  the 
donor  energy  levels  induced  by  the  random  distribution  of  electric  fields 
at  the  sites  of  donors  in  the  crystal.  These  fields  originate  from  the 
spatial  distribution  of  charged  acceptors  and  ionized  donors.  D.M.  Larsen 
has  predicted  line  shapes  for  the  ls-2p  transitions  on  the  basis  of  such  a 
model  which  agree  well  with  line  shapes  observed  in  high  purity  GaAs  samples.^ 
In  very  pure  uncompensated  samples  when  the  field  has  reached  50  kilogauss. 
the  ls-2p  (m=  -1)  line  width  has  shrunk  from  a  zero  field  value  of  about 
5  cm  ^  to  as  little  as  0.5  cm  ^  and  peaks  corresponding  to  transitions  to 
other  hydrogenic  excited  states  are  clearly  resolved.  Figure  4  shows  a 
high  resolution  spectrum  of  a  high  purity  sample  at  52.5  kilogauss.  Notice 
that  both  the  ls-2p  (m=  -1)  and  the  ls-2p  (m=Q)  transitions  have  developed 
a  fine  structure  which  was  not  apparent  at  zero  field.  This  fine  structure 
is  the  result  of  the  central  cell  shift  of  the  donor  ground  state  energy, 
which  is,  in  general,  different  for  different  impurity  species.  The  energy 
of  the  2p  state  is  almost  totally  unaffected  by  thu  central  cell  potential 
since  the  p  wavefunctions  are  zero  at  r*0.  As  a  result,  the  is-2p  transition 
energy  is  impurity  dependent  primarily  through  the  Is  central  cell  correction. 

The  ls-2p  (m*-l)  transitions  are  of  particular  interest  since  they  have 
the  largest  amplitude  and  are  the  narrowest  in  .e  50  kilogauss  spectrum. 

At  this  high  field, three  conditions  have  conspired  to  facilitate  impurity 


RELATIVE 


RELATIVE  PHOTO  CONDUCT!  VE  RESPONSE 


r 


Figure  3.  Photothermal  ionization  spectra  at  several  magnetic  fields 
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Figure  A.  High  resolution  phototherraal  ionization  spectra  in  GaAs. 


identification  using  the  ls-2p  (m^-l)  transition.  The  ls-2p  (m— 1)  lines 

have  been  Zeeman  split  away  in  energy  from  the  other  broader  ls-2p  lines. 

The  electric  quadrunolc  moment  of  the  2p  (m— 1)  state  has  been  reduced 

nearly  to  zero  and  so  the  main  contribution  to  Stark  broadening  of  the 

ls-2p  (m*l)  transition  has  been  removed.  Finally, the  compression  of  the 

Is  wavefunction  by  the  magnetic  field  has  increased  the  overlap  with  the 

central  cell  and  so  has  increased  the  central  cell  shifts  of  the  Is  state. 

As  a  result  in  very  high  purity  GaAs, it  is  possible  to  resolve  peaks  in  this 

transition  corresponding  to  many  different  donor  species.  The  relative 

position  of  a  peak  within  the  ls-2p  (m=-l)  envelope  has  been  used  as  a 

signature  of  a  given  donor  species,  and  many  of  the  donors  in  GaAs  have 

been  tentatively  identified  in  this  way.  The  identifications  due  to  C.M. 

3 

Wolfe,  G.E.  Stillman  and  D.M.  Korn  are  shown  in  Fig. 5. 

A  block  diagram  of  the  apparatus  used  for  the  far  infrared  photo¬ 
conductivity  measurements  appears  in  Fig. 6(a).  It  consists  of  a  spectrometer, 
a  cryostat,  and  the  LSI-11  minicomputer  based  data  acquisition  system.  The 
cryostat  Includes  a  superconducting  magnet  which  can  expose  the  sample  to  a 
0.1%  uniform  field  of  up  to  67  kilogauss.  Because  the  spectrometer,  coupling 
optics,  and  cryostat  are  sensitive  to  vibrations,  they  have  been  mounted  on 
an  optical  table  which  in  turn  rests  on  a  large  concrete  pillar  sunk  in  the 
ground  (separate  from  the  building's  floor  vibrations). 

The  spectrometer  (see  Fig. 6(b))  is  a  Beckman  FS720  Michelson  interferometer 
which  is  used  as  a  high  resolution  (±0.06  cm  ^)  tunable  (5  cm  ^  to  1000  cm  ^) 
far  infrared  light  source  for  illuminating  the  semiconductor  sample.  Its 
operation  can  be  understood  in  the  following  way.  The  beamsplitter  divides 
the  collimated  light  from  the  source  into  two  beams.  One  beam  travels  a 
fixed  path  length  while  the  other  travels  a  path  length  which  is  variable 
by  means  of  the  spectrometer's  movable  mirror.  At  the  beamsplitter  the  two 
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Figure  5a.  ls-2p  (m«  -1)  peak  positions  at  50  kilogauss 
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Stillman  et  al.  identification 


*  For  each  run  the  peaks  have  been  shifted  so  that  the  Si  peak  agrees  with  Stillman's 
identification  to  account  for  slight  variations  in  magnetic  field  from  run  to  run. 
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Figure  ti  (b) 

Figure  6.  Apparatus  block  diagram. 


beams  recombine  and  are  directed  to  the  sample  by  means  of  light  pipes. 
When  the  path  difference  is  zero,  the  two  beams  have  the  same  phase  for  all 
wavelengths  and  the  recombined  beam's  intensity  is  a  maximum.  For  a  mono¬ 


chromatic  source  at  wave  number  v  and  beams  of  equal  intensity  S_,  the  combined 
beam  intensity  as  a  function  of  path  difference  x  is  given  by 

(4)  S(x)  =  S-  (1+cos(2ttvx)) 

The  photoconductivity  signal  at  the  sample  is  then  given  by 

(5)  Iv  (x)  =  R(u)  [T-S-(l+cos (2ttvx)  )  ] 

where  T_  is  the  total  transmission  of  the  spectrometer  optics  (filters,  beam¬ 
splitter  efficiency,  etc.)  and  R(v)  is  the  object  of  interest,  the  responsivity 
of  the  sample  at  wavenumber  u.  Since  the  source  emits  a  continuous  spectral 
distribution  S(\5)  the  total  photoconductivity  signal  is  given  by  an  integral 
over  all  wavenumbers. 

00 

I(x)  =  /  R(v)S(v)T(v)  [1+cos(2ttvx)  ]dv 
o 

The  Fourier  cosine  transform  of  this  expression  gives: 

00 

(6)  R(S)S(")T(v)  =  2  /  [I(x)-I(0)/2]e'i2^Xdx 

— GO 
00 

=  4  /  [I(x)-I(0) /2]cos(2iTvx)dx 
o 

The  interferogram  I(x)  is  the  quantity  measured  experimentally  by  recording 
photoconductivity  signal  as  a  function  of  mirror  position.  The  product 
S(v)T(v)  can  be  established  using  (6)  by  a  measurement  taken  with  a  detector 
of  known  or  flat  responsivity.  The  responsivity  of  the  semiconductor  sample 
R(v)  then  is  essentially  the  Fourier  transform  of  the  recorded  interferogram 


To  record  the  interferogram,  the  biased  sample  is  illuminated  with 
light  from  the  spectrometer  which  is  chopped  at  150  Hz.  The  ac  component 


of  the  voltage  across  the  sample  is  proportional  to  the  difference  in 

sample  conductivity  with  the  light  on  from  that  with  the  light  off. 

This  signal  is  amplified  and  fed  to  a  lock-in  amplifier  referenced  with  a 

signal  from  the  spectrometer's  chopper.  Each  time  the  LSI-11  interface 

receives  a  pulse  from  the  Moir6  grating  electronics  in  the  spectrometer 

(indicating  that  a  specified  Ax  in  mirror  position  has  been  traversed) 

the  LSI-11  records  the  lock-in  output  through  an  A/D  converter  as  one  data 

point  in  the  interferogram.  The  spectrometer  mirror  is  in  continuous  motion 

at  speed  v  and  the  lock-in  amplifier's  time  constant  is  chosen  to  be  comparable 

to  the  time  between  digitize  pulses  Ax/v.  The  recorded  lock-in  output  then 

is  an  average  of  the  photoconductivity  amplitude  over  the  last  Ax  of  mirror 

travel.  Evidently  from  equation  (A),  Ax  is  to  be  chosen  to  be  at  least 

as  small  as  1/(2  v  )  where  v  is  the  largest  wavenumber  of  interest.  When 
max  max  ° 

a  certain  predetermined  number  of  data  points  have  been  recorded  (typically 
4096),  the  completed  interferogram  is  stored  on  the  LSI-11 's  floppy  disc. 

The  interferogram  can  be  quickly  displayed  on  a  graphics  terminal,  or,  if 
hard  copy  is  desired,  it  can  be  permanently  recorded  on  an  x-y  plotter. 

An  example  of  such  an  interferogram  plot  is  shown  in  Fig. 7. 

After  it  is  recorded,  mathematical  operations  on  the  interferogram  data 
can  be  performed  by  the  LSI-11  leading  to  a  spectrum  of  photoconductive 
response  versus  wavenumber.  These  operations  may  include  removal  of  noise 
spikes,  slope  or  D.C.  shift  in  the  interferogram,  apodization,  zero  filling, 


etc.,  but  the  essential  operation  is  a  discrete  Fourier  transform  (DFT) 
of  the  interferogram. 


transform  spectrum  from  the  apparatus 


The  photoresponse  of  the  sample  R ( \5)  is  related  to  the  interferogram 


I(x)  through  equation  (6). 

00 

(7)  R(v)S(v)T(v)  =  4  /  [I(x)-I(0) /2]cos(2irvx)dx 


=  4  Re{/  [I(x)-I(0)/2]e-i27rvxdx} 


Since  all  calculations  are  to  be  performed  by  a  digital  computer  (the  LSI-11 


in  our  case),  I(x)  is  sampled  at  discrete  values  x^  over  a  finite  interval 


(x^  =  nD/(N-l)  n  =  0,  1,  .  .  .  N-l)  and  a  discrete  approximation  to  the 


above  integral , 


^  _  -I  VrjY 

(8)  R(  v) S (v) T(v)  -  4  Re{T  [I(x  )-I(0)/2]e  ^  n  Ax  } 

u  n  n 

n=o 


N-l 


=  4 (D/ (N-l) ) Re{ \  [I(xn)-I(0)/2]e"i2ffvxn> 

n=o 


must  be  used. 


The  quantity  R(v)S(v)T(v)  is  calculated  only  at  N  discrete  values 


—  t2tt  /N 

(v^=k/D  k=0,  1,  .  .  . ,N-1).  With  the  definition  e  HW  ,  ecluat:*-on  (8) 


can  be  written  as 


N-l 


(9) 


R(\)3(vk)T(vk)  =  4  (D/(N-l)Re{£  [I(xn)-I(0) /2]wJ 


nk. 


n=o 


If  we  define  [I(x  )-I(0)/2]  -  a(n)  and 
n 


R(vk)f (Vk)T(vk)/(4D/(N-l))  H  Re[A(k) 


then  (9)  can  be  written 


N-l 


(10) 


A(k)  =  l  a(n) 


,+nk 
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Thus  A(k),  which  is  proportional  to  the  photoresponse  of  the  sample,  is  the 

DFT  of  the  interferogram  a(n) .  The  inverse  transform  is  just 

N-l 

(ID  a(n)  -  |  l  A(k)W~nk  • 

k=o  u 

It  has  been  tacitly  assumed  in  all  of  the  above  discussion  that  the  point 

where  the  path  difference  in  the  spectrometer  is  zero  is  x  =0  which  corres- 

n 

ponds  to  the  zeroth  sampled  point  in  the  interferogram.  This  is  not  the 
case  experimentally.  The  mirror  position  of  zero  path  difference  (ZPD)  will 
lie  at  some  xZpQ: ®-xzpDiD  and  in  general  will  lie  between  sampled  inter¬ 
ferogram  points  (See  Fig. 8).  Thus  the  analysis  should  proceed  from  equation 
(4)  being  changed  to  read: 


S(x)  =  S_  (l+cos2^v(x-xzpD)) . 


For  a  continuous  Fourier  transform  the  fact  that  xZpD^0  presents  no  problem 


since  by  a  change  of  variables 


/  I(x-x7pn)e-i27^xdx  -  e‘i2™xZPD  /  W*'^*** 


(12) 


ZPD 


We  can  then  correct  for  the  fact  that  zero  path  difference  was  not  at  x=0 
simply  by  knowing  *ZpD  and  multiplying  by  an  exponential  factor.  However, 
in  sampling  the  interferogram  discretely,  if  a  measurement  is  not  taken  at 
XZPD  *'as  :''S  8enerally  c^e  case),  serious  distortions  can  appear  in  the 


computed  DFT  spectrum.  Such  distortion  is  representative  of  one  class  of 
phase  error,  and  there  are  a  number  of  phase  correction  schemes  designed  to 
eliminate  this  problem. 


Shown  is  an  example  of  a  hypothetical  interf erogram  in  which  no  sample 
point  occurs  at  zero  path  difference.  Note  the  resulting  asymmetry  in 
the  digitized  interf erogram  below. 


Figure  8.  Graphic  representation  of  a  course  of  error. 


Two  such  phase  correction  schemes  have  been  applied  by  us  to  both 
experimentally  obtained  and  synthetic  interf erograms .  Their  effects  on 
spectral  line  widths  and  line  shapes  have  been  compared  for  a  variety 
of  different  synthetic  interf erograms . 

The  first  phase  correction  method  is  extremely  simple  and  can  be  under¬ 
stood  by  examining  equation  (12).  Since  the  interf erogram  is  symmetric 
about  xzpD  the  Fourier  transform  with  respect  to  an  origin  at  xzpD  is  real 
and  equal  to  the  modulus  of  either  side  of  equation  (12).  Thus  the  first 
means  of  phase  correction  consists  of  taking  the  full  complex  DFT  of 
[I(xn)-I(0) /2] ]  and  then  taking  the  modulus  of  the  result. 

Another  method  of  phase  correction  we  have  studied  is  similar  to  that 

g 

described  by  Mertz.  In  this  scheme  the  full  complex  DFT  of  the  severely 
apodized  interf erogram  is  computed  and  the  phase  of  this  specturm  is 
recorded.  Severe  apodization  is  accomplished  by  multiplying  the  inter- 
ferogram  [ (I(xn)-I(0) /2]  by  a  function  a g(x)  of  the  form  (See  Fig. 9(a).) 

I x-x  I 

(13)  a  (x)  =  1-  - for  |xj<xA 

s  x.  ' "  A 

A 

=  0  for  | x  j  >x . 

where  x&<<D  and  xmax  is  the  position  of  the  interf  erogram  maximum  (which 

is  presumably  within  D/(N-1)  of  x^pp).  The  phase  $g  constitutes  an  average 

or  reference  phase  indicative  of  the  position  of  the  true  zero  path  difference 

XZPD‘  Then  the  full  complex  DFT  of  the  normally  apodized  interf erogram  is 

computed.  Normal  apodization  is  accomplished  by  multiplying  the  interferogram 

by  a  (x) .  (See  Fig.  9(b).)  This  functional  form  of  a  (x)  is  chosen  so  that  the 
11  n 

redundant  information  on  either  side  of  x  sx__„  is  not  double  counted. 

max  ZPD 

Finally,  the  phase  corrected  spectrum  Acorrected 00  is  computed  by 
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Apodization  Functions 


Figure  9  (b) 

Figure  9.  Apodization  functions. 


(14) 


Corrected  (k)  =  Re^(k)e“i$s}  . 

It  can  be  shown  that  if  a  (x)  were  chosen  to  be  a  symmetric  box  about  x 

s  max 

then  the  method  described  above  is  exactly  equivalent  to  the  method  described 
by  Mertz. 

For  the  study  of  impurity- impurity  interactions  through  line  shapes  in 
the  photothermal  ionization  spectra,  it  is  essential  to  ascertain  and  minimize 
any  effect  of  the  digital  Fourier  transformation  technique  on  those  line  shapes. 
In  order  to  compare  the  two  phase  correction  schemes  described  above,  and 
in  particular  their  effect  on  spectral  line  shape,  a  program  was  written  to 
generate  synthetic  interferograms .  These  are  digitally  sampled  finite  path 
difference  subsets  of  a  continuous  function  whose  Fourier  integral  transform 
is  a  wedge  of  the  form 


(15) 

f(v) 

■  0 

v<vx 

f(?) 

*  +A(v-v1) /(v2-v^) 

f(v) 

*  -A(v-v3)/(v3-v2) 

V2iV^3 

f(v) 

-  0 

v.<v 

The  form  of  this  function  was  chosen  to  simulate  the  characteristic  line  shape 
observed  in  the  ls-2p  (m=-l)  transitions  in  GaAs.  The  synthetic  interferograms 
were  generated  using  various  sampling  parameters  typical  of  those  used  experi¬ 
mentally.  Also,  and  most  important  from  the  standpoint  of  comparing  phase 
correction  methods,  interferograms  were  generated  with  various  amounts  of 
phase  error.  That  is,  the  distance  from  the  continuous  interferogram's 
maximum  to  the  nearest  sampled  point  was  varied.  (See  Fig.  8.)  The  discrete 


Fourier  transforms  of  these  synthetic  interf erograms  were  then  computed 
and  the  resulting  spectra  compared. 

In  an  interferogram  from  a  well  aligned  spectrometer,  the  data  are 
symmetric  about  xzpD*  C^e  ^nf°mation  *-s  contained  in  one  or  the 

other  halves  of  the  interferogram  and  so  recording  both  halves  is  a  waste 
of  time  and  computer  memory.  Also,  the  resolution  Av  (cm  of  the 
computed  spectrum  is  given  by 

(16)  Av  =  i. 

where  D  (cm)  is  the  maximum  optical  path  difference  achieved  by  the  spectro¬ 
meter's  movable  mirror  during  the  recording  of  the  interferogram.  In  high 
resolution  spectra,  to  minimize  Av,  D  is  chosen  as  large  as  possible.  The 
mirror  drive  in  the  spectrometer  can  move  the  mirror  a  maximum  of  10  cm 
which  corresponds  to  an  optical  path  difference  change  of  20  cm.  To  improve 
resolution,  the  spectrometer  has  been  modified  in  such  a  way  as  to  move 
from  a  mirror  position  near  the  middle  of  the  mirror  drive  travel  to  a 
position  near  one  end.  This  allows  the  recording  of  one-sided  interf erograms 
with  a  D  of  nearly  20  cm  and  with  a  corresponding  resolution  of  about  0.05  cm  . 
With  these  one-sided  inter f erograms ,  proper  phase  correction  is  a  necessity. 

Our  synthetic  interferograms  have  been  constructed  to  simulate  the  experi¬ 
mentally  recorded  interferograms,  having  xZpD=0.A8  cm  and  D=18.0  cm.  The 
number  of  sample  points  N  and  sample  interval  Ax  were  N=3850  and  Ax=48  microns. 

Table  III  summarizes  the  sampling  parameters  and  phase  error  information 
and  figures  10  through  15  show  some  of  the  resulting  soectra.  The  labels 
"modulus  spectrum"  and  "phase  referenced  spectrum" indicate  that,  respectively, 
the  first  and  second  of  the  above  described  phase  correction  methods  were 
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TABLE  III.  SYNTHETIC  INTERFEROGRAM  PARAMETERS 


Figure  10.  Modulus  spectrum  of  synthe 


Phase  referenced  spectrum  of  synthetic  interf eropram 


spectrum  of  synthetic  inter ferogram 
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employed  in  computing  that  spectrum.  Notice  that  in  all  these  spectra,  as 
in  the  experimental  spectra,  the  line  widths  are  extremely  narrow,  being 
from  .014  to  .007  of  the  wavenumber  at  which  they  occur. 

Figures  10  and  11  are  respectivelv  the  modulus  and  phase  referenced 
transforms  of  inter ferogram  0  which  has  v^=35.0  cm  ^ ,  ~^=36.0  cm  \  and 
u^=36.5  cm  A  perfect  undistorted  transform  of  this  interferogram  would 
be  a  triangular  wedge  with  vertices  at  these  wavenumbers,  and  both  spectra 
approximate  this.  The  phase  error  for  this  interferogram  was  0.5  Ax  which 
is  considered  a  worst  case. 

Notice  that  the  modulus  spectrum  is  the  poorer  of  the  two  approximations 
to  the  true  line  shape,  having  too  large  a  line  width,  distortion  of  the  ideally 
straight  sides  of  the  triangular  wedge,  and  broad  side  lobes  which  extend 
out  several  wavenumbers  from  the  main  peak.  These  side  lobes  are  not  the  ones 
usually  associated  with  unapodized  spectra,  having  a  period  which  is  far 
too  large  to  be  accounted  for  in  this  way. 

The  phase  referenced  spectrum  is  a  considerably  closer  approximation 
to  the  ideal  line  shape  but  it  too  is  slightly  distorted.  The  side  lobes 
apparent  in  the  modulus  spectrum  are  absent  here,  but  the  distortion  in  the 
main  body  of  the  peak  is  similar  in  character  though  much  reduced.  The 
line  width  in  the  phase  referenced  spectrum  is  correct  and  the  only  distortion 
away  from  the  main  peak  is  a  slight  upward  curvature  of  the  baseline  as  the 
peak  is  approached  from  either  side. 

The  line  width  of  interferogram  0  is  considerably  larger  than  one 
encounters  in  the  purest  C.aAs  available.  In  Fis.  22  for  example,  the  spectrum 
of  B61-11  has  line  widths  of  order  0.1  cm  \  With  this  in  mind,  interferogram  7 
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was  prepared  with  v^=35.75,  V2=36.00,V2=36.00,  and  a  phase  error  of  0.5  Ax. 

The  modulus  and  phase  referenced  spectrum  of  this  interferogram  are  shown 
respectively  in  figures  12  and  13. 

Again  the  modulus  spectrum  represents  the  poorer  approximation  to  the 
true  line  shape.  Its  line  width  is  too  large  and  the  ideal  precipitous  drop 
in  amplitude  at  high  wavenumber  is  smeared  into  a  gradual  decline.  The 
sidelobes,  although  not  apparent  in  the  figure,  were  present  with  amplitude 
of  about  0.02  that  of  the  main  peak.  There  is  a  pronounced  peaking  at  the 
apex  of  the  line  and  a  gross  distortion  of  the  ideal  right  triangular 
peak  shape. 

The  phase  referenced  spectrum  is  a  much  closer  approximation  to  the 
true  line  shape.  The  most  obvious  deviations  are  the  Gibbs  oscillations 
around  the  high  wavenumber  edge  of  the  peak  which  were  mentioned  above. 

The  ideal  precipitous  drop  in  amplitude  at  ~2  is  again  smeared,  but  this  time 
only  over  a  wavenumber  range  equal  to  the  resolution  Av=.05  cm  .  The  low 
energy  side  of  the  peak  is  slightly  concave  from  above  and  contains  the  Gibbs 
oscillations,  but  it  is  a  fairly  good  approximation  to  the  true  line  shape. 
Interferogram  7  was  prepared  identically  to  interferogram  6  except  that  x^p^ 
was  placed  in  the  middle  of  interferogram  7.  The  modulus  spectrum  of  7  is 
shown  in  Fig.  14.  Because  the  distance  from  t0  extremes  °f  the  inter¬ 

ferogram  has  been  halved,  the  resolution  of  the  spectrum  is  likewise  halved, 
but  it  is  interesting  to  note  that  the  line  shape  and  line  width  are  nearly  as 
good  as  those  of  the  phase  referenced  spectrum  of  interferogram  6. 

Finally , inter ferograms  10,  11,  12  and  13  were  prepared  to  investigate 
the  size  and  nature  of  the  effect  of  phase  error  on  lineshape.  All  four 
inter f erograms  are  identical  having  v^=35.0,  V2~36.0,  and  v^=36.0  except  for 
their  phase  errors  which  in  units  of  the  sampling  interval  Ax  are  0.05,  0.1, 
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0.3  and  0.5  respectively.  Figure  15  is  a  Dlot  of  the  nhase  referenced  spectra 

of  these  interf erograms  which  are  labeled  respectively,  1,  2,  3  and  4.  As 

the  phase  error  is  increased,  that  is, as  the  closest  sample  point  to 

is  moved  progressively  further  away  from  x  ,  the  distortion  of  the  true 

ZPD 

line  shape  becomes  worse  until  at  a  phase  error  of  0.5  it  is  at  a  maximum. 

The  maximum  deviation  from  the  ideal  line  shape  is  about  10%  of  the  peak 
amplitude  and  this  might  well  be  significant  if  one  were  examining  line  shapes 
to  study  impurity  interactions.  We  are  currently  investigating  possible 
means  for  reducing  this  distortion  further. 

EXPERIMENTAL  WORK 

Several  samples  of  high  purity  GaAs  grown  by  various  industrial  and 
government  laboratories  have  been  measured  using  the  photothermal  ioniza¬ 
tion  technique.  These  samples  are  representative  of  the  best  quality 
material  which  can  presently  be  grown  by  various  growth  techniques.  Hall 
effect  measurements  at  77  K  in  a  field  of  6  kilogauss  have  determined  the 
mobility  and  donor  and  acceptor  impurity  concentrations  (N^  and 

N , ,  respectively)  in  all  of  these  samples.  (See  Table  IV.)  3v  examining 
A 

the  ls-2p(m=-l)  peak  positions  in  the  50  kilogauss  spectra  of  these  samples 

3 

and  comparing  with  the  identit ications  ot  Wolfe,  Stillman,  and  Korn  ,  we 
hc.ve  established  the  residual  donor  species  present  in  these  samples.  The 
correlation  between  residual  donors  present  and  growth  technique  for  samples 
grown  in  different  laboratories  suggests  that  the  residual  donors  we  have 
identified  may  be  characteristic  of  the  associated  growth  technique. 

MBE  SAMPLES 

Three  samples  grown  by  molecular  beam  epitaxy  (MBE)  were  measured.  The 

first  of  these,  MBE  A-147,was  grown  by  A.Y.  Cho  and  H.  Morko?  at  Bell  Labs. 

3  5 


The  remaining  samples  MBE-91  and  MBE-92  were  grown  by  A.R.  Calawa  at  MIT 
Lincoln  Labs.  Photothermal  ionization  spectra  for  these  samples  appear  in 
Figures  16,  17  and  18. 

The  spectrum  of  MBE  A-147  shows  three  distinct  peaks  which  can  be  associated 
with  known  simple  donor  species.  The  first  and  largest  of  these  is  either  Se 
or  Sn  (the  ls-2p  (m=-l)  transition  energies  for  these  two  species  being  very 
similar) .  A  second  peak  of  much  smaller  line  width  occurs  at  higher  wavenumber 
and  is  almost  certainly  Si.  Finally,  at  a  still  higher  wavenumber,  a  small  peak 
attributable  to  C  can  be  discerned  above  the  noise.  The  relative  amplitudes 
of  these  peaks  give  a  rough  indication  of  the  relative  concentration  of  the 
associated  donor  species.  (The  integrated  area  under  each  peak  would  be 
a  better  measure.)  These  peak  amplitudes  are  in  a  ratio  of  (SeSn) :Si:C: : 
20.0:8.5:22.0.  It  is  interesting  to  note  the  marked  difference  in  line  width 
between  the  (SeSn)  peak  and  the  Si  peak  since  in  the  usual  picture  of  a 
uniform  distribution  of  donor  and  acceptor  impurities  the  ls-2p  (m=-l) 

Stark  broadening  is  expected  tc  be  the  same  for  all  donor  species.  This 
suggests  the  possibilities  of  either  an  additional  source  of  broadening  or 
a  donor  species  dependent  non-uniform  distribution  of  charged  impurities 
about  the  donor  site  or  an  otherwise  inhomogeneous  sample  7. 

The  spectra  of  MBE-91  and  MBE-92  show  the  same  residual  impurities  to 

be  present  but  in  rather  different  relative  concentrations.  Again  as  in 

MBE  A-147, the  (SeSn)  and  Si  peaks  are  prominent  but,  in  addition, Pb  and 

possibly  X^  (seen  as  a  shoulder  on  the  (SeSn)  peak)  are  present.  The  X^ 

peak  is  seen  in  VPE  material  as  a  residual  donor  but  to  our  knowledge  it 

has  never  been  seen  in  LPE  material.  It  has  been  suggested  that  X^  may  be 

12 

a  Ga  vacancy- impurity  complex  but  this  identification  must  be  regarded 
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as  speculative.  The  line  widths  of  all  the  peaks  in  both  spectra  are 
comparable  and  the  relative  peak  amplitudes  are  for  MBE-91:  Pb:X^: (SeSn) : 
Si: :58. 8:77. 0:78. 5:16. 6,  and  for  MBE-92:Pb:X1: (SeSn) :S1: :69. 5:33. 0:43. 2: 
33.4. 

The  presence  of  Pb  and  X^  in  MBE-91  and  MBE-92  turns  attention  once 
again  to  the  abnormally  broad  peak  labeled  (SeSn)  in  the  spectrum  of  MBE 
A-147.  It  is  conceivable  that  the  width  of  this  peak  may  be  due  to  the 
presence  of  unresolved  Pb  and/or  X^  peaks  at  lower  energies,  and  in  fact 
small  shoulders  on  the  (SeSn)  peaks 's  low  energy  side  can  be  seen  at 
wavenumbers  associated  with  these  species  but  these  are  too  small  to  be 
considered  definitive. 

In  summary,  all  the  MBE  grown  material  has  peaks  associated  with 
(SeSn)  and  Si.  Both  Pb  and  X^  occur  in  the  MIT  samples  and  possibly 
in  the  Bell  Labs  sample,  while  C  is  present  only  in  low  concentration  in 
the  Bell  Labs  sample  and  absent  entirely  in  the  MIT  samples. 

MO-CVD  SAMPLES 

P.D.  Dapkus  of  Rockwell  International  has  provided  us  with  three 

GaAs  samples  MO-CVD  Seq.  29,  Seq.  129,  and  Seq.  132  grown  by  the  metal 

organic  chemical  vapor  deposition  (MO-CVD)  growth  technique.  Photo- 

thermal  ionization  spectra  for  these  samples  aopear  in  Figures  19,  20,  and 

21.  In  all  three  spectra,  the  dominant  neak  corresponds  to  that  identified 

3 

by  Wolfe,  Stillman  and  Korn  as  C.  That  C  should  be  incorporated  into 
the  MO-CVD  grown  material  is  not  surprising  in  view  of  the  organic  carrier 
molecules  present  during  growth.  The  dominance  of  this  peak  in  all  three 
spectra  can  be  seen  as  a  confirmation  of  the  identification  of  C  by 
Wolfe  et  al.,  as  well  as  a  determination  of  the  primary  residual  donor 


species  in  MO-CVD  material. 
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-CVD  Seq.129 


38.00 


It  should  be  noted  here  that  M.N.  Asfar,  K.S.  Button,  and  G.L.  McCoy, 
using  neutron  transmutation  doping,  have  observed  a  peak  at  approximately  the 
position  (relative  to  the  other  peaks  in  the  ls-2p  (m=-l)  envelope)  of  the 
peak  we  have  labeled  as  C.  Asfar  et  al.  have  identified  their  peak  as  Ge  in 
accordance  with  the  expected  capture  of  thermal  neutrons  by  Ga  nuclei  in  the 
GaAs  lattice. 

Ga  +  n  -*e  +Ge 

As  +  n  ->-e  +Se 

The  corresponding  neutron  capture  reaction  by  As  results  in  transmutation  to 

Se  whose  photothermal  ionization  peak  position  has  been  identified  by  C.M. 

3 

Wolfe,  G.E.  Stillman,  and  D.M.  Korn.  Asfar  et  al.  also  observed  a  peak  at  a 
position  in  agreement  with  this  identification  in  their  transmutation  doped 
material. 

In  the  text  of  their  paper,  Asfar  et  al.  suggest  that  the  position  of 

the  peak  they  identify  with  Ge  should  not  be  associated  with  C.  There  is  no 

reason  to  suspect  a  priori  that  a  peak  position  should  be  uniquely  held  by  a 

single  donor  species.  This  is  especially  true  in  view  of  the  large  number  of 

donor  species  and  the  narrow  range  of  wavenumbers  (compared  to  a  peak  line 

width)  over  which  the  ls-2p  (m=-l)  transitions  occur.  Moreover,  it  is 

difficult  to  imagine  how  the  dominant  peak  in  the  spectra  of  all  of  our  MO-CVD 

samples  could  be  attributed  to  Ge  and  not  C.  In  addition,  the  Ge  identification 

3 

by  Asfar  et  al.  contradicts  identif icat ions  by  C,  Wolfe,  et  al.,  R.  Cooke 
et  al.,^  and  M.  Ozeki  et  al.^  All  three  of  these  groups  have  identified  a 
peak  at  lower  wavenumber  (about  O.b  cm  1  above  the  Si  peak  at  50  kilogauss) 
in  Ge  back  doped  samples.  Perhaps  Ge  introduces  more  than  one  shallow  donor 
level  in  GaAs.  If  C  and  Ge  do,  in  fact,  introduce  peaks  at  or  near  the  same 
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wavenumber  within  the  ls-2p  (m=-l)  envelope, this  may  hinder  or  even  preclude 
discrimination  between  these  donor  species  in  GaAs. 

In  the  spectrum  of  MO-CVD  Seq.  29,  a  small  peak  at  a  wavenumber  corresponding 
to  Si  is  resolved,  while  in  the  spectra  of  Seq.  129  and  Seq.  132  the  Si  peak 
has  more  than  half  the  amplitude  of  the  dominant  C  peak.  Also  in  the  spectrum 
of  Seq.  132  (which  was  grown  at  600°C)  a  large  peak  at  the  position  of  X^ 
is  seen,  while  in  the  spectrum  of  Seq.  129  (which  was  grown  at  575°C  under 
conditions  otherwise  identical  to  those  of  Seq.  132)  the  X^  peak  is  absent. 
Because  of  the  possible  association  of  X^  with  a  Ga  vacancy- impurity  complex, 
this  may  be  an  indication  of  increased  vacancy  formation  at  the  higher  growth 
temperature. 

The  relative  peak  amplitudes  in  these  samples'  spectra  are  Si:C:: 

1.8:6. 8  in  Seq.  29,  Si:C: :6.0:1.04  (at  25  kilogauss)  in  Seq.  129,  and 
X1:Si:C::3. 54:4. 40:7. 06  in  Seq.  132. 

VPE  SAMPLES 

Three  samples  grown  by  vapor  phase  epitaxi  (VPE)  were  measured.  The 
first  of  these  B61-11  was  grown  by  C.M.  Wolfe,  while  he  was  at  MIT  Lincoln 
Laboratory  using  an  AsCl^  growth  technique.  The  other  two  samples  GaAs- 
763  and  GaAs  775  were  grown  by  J.K.  Kennedy  at  Hanscom  Air  Force  Base  using 
the  AsH^  hydride  technique.  The  spectra  of  all  three  samoles  have  neaks 
due  to  X^,  Si,  and  C.  (See  Figs.  22,  23,  and  24.) 

The  spectrum  of  B61-11  shows  the  extremelv  small  line  widths  associated 

with  such  high  purity  and  lightly  compensated  material.  Donor  and  acceptor 

concentrations  in  this  sample  were  determined  from  77  K  Hall  measurements 

cm  ^  and  N  =9.44xlo'"‘  cm  ^  with  t__  =  200,993  cm^/Vsec. 
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This  is  the  best  quality  GaAs  sample  we  have  measured,  and  because  of  the 
narrow  peaks  in  its  photoconductivity  spectrum,  it  has  been  used  as  a 
standard  for  identifying  peaks  in  the  spectra  of  other  samples.  The 
oscillations  seen  in  the  spectra  of  B61-I1  are  mostly  due  to  the  convolu¬ 
tion  of  the  extremely  narrow  lines  with  the  instrument  function  of  the 
spectrometer.  Alternatively, they  can  be  thought  of  as  Gibbs  phenomenon 
resulting  from  the  finite  path  length  over  which  the  interf erogram  must 
be  measured.  These  oscillations  decrease  in  wavelength  and  amplitude  as 
the  path  length  over  which  the  interf erogram  is  measured  is  increased. 

The  spectrum  of  B61-11  is  typical  of  that  of  the  best  undoped  samples 
grown  in  C.M.  Wolfes'  AsCl^  system  except  that  the  relative  amplitude  of  the 
X^  peak  is  a  bit  larger  and  of  the  C  peak  is  a  bit  smaller  with  respect  to 
the  Si  peak  than  has  been  usually  seen.  The  ratios  of  peak  amplitudes 
in  the  spectrum  of  B61-11  are  X^:Si:C: : 8 3 . 1:78. 6:17. 8, while  in  eight  other 
high  purity  undoped  GaAs  samples  grown  in  this  system, spectra  taken  by 
G.E.  Stillman  at  MIT  Lincoln  Laboratory  show  average  relative  peak  amplitude 
of  X1:Si:C::0.22:1.0:0.76. 

In  the  spectra  of  both  the  samples  GaAs-763  and  GaAs-775,  only  peaks 
due  to  X^,  Si  and  C  are  seen.  The  relative  peak  amplitudes  in  GaAs-763  are 
X1:Si:C::62. 3:312:11. 3  and  in  GaAs-775  are  X  :Si:C: :41. 6 : 94 . 3 : 13 . 5 .  Notice 
that  in  both  these  spectra  the  Si  peak  dominated  followed  in  amplitude  by 
the  X^  and  then  the  C  peak.  The  peaks  in  the  spectrum  of  GaAs-763  are 
considerably  narrower  than  those  of  GaAs-775  as  would  be  expected  from  the 
factor  of  two  difference  in  acceptor  concentration  which  would  cause 
difference  in  intensity  of  random  electric  fields  responsible  for  Stark 
broadening. 


LPE  SAMPLES 


L.F.  Eastman  of  Cornell  University  has  provided  us  with  two  GaAs 
samples  SM-25A  and  SM-30  grown  by  the  liquid  phase  epitaxy  (LPE)  growth 
technique.  Another  sample  LPE-251(2)  grown  by  J.  Abrokwa  at  Honeywell 
using  this  technique  was  also  measured.  Photothermal  ionization  spectra 
for  these  samples  appear  in  Figures  25,  26  and  27. 

For  the  spectra  of  both  SM-25A  and  SM-30  peaks  due  to  either  Se  or 
Sn,  and  to  Si  are  observed  with  roughly  equal  amplitudes.  As  was  the  case  in 
the  two  hydride  VPE  samples,  the  LPE  sample  SM-25A,  having  lower  acceptor 
concentration,  has  less  Stark  broadening  in  the  peaks  of  the  photoconductivity 
spectrum.  (See  Table  III  for  the  values  of  and  line  widths.) 

The  sample  LPE-251(2)  is  interesting  in  that  its  spectrum  shows  only 
one  donor  species  is  present.  The  peak  in  this  spectrum  is  that  of  Si,  and 
although  there  is  some  structure  on  the  low  energy  side  of  this  peak  at 
wavenumbers  close  to  those  associated  with  X^  and  (SeSn),  such  structure 
may  well  be  noise  in  the  spectrum.  It  would  be  surprising  to  find  an  X^ 
peak  because  to  our  knowledge  this  has  never  been  seen  in  LPE  grown 
material.  Another  interesting  aspect  of  this  spectrum  is  the  narrowness 
of  the  Si  peak.  This  is  incommensurate  with  the  rather  low  mobility  and 
high  acceptor  concentration  deduced  for  this  sample  from  our  77  K  Hall 
measurement . 
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